% 计算切比雪夫多项式基函数
% 从下标0开始基函数
% 输入参数：
% N: 基函数的个数
% t: 时间变量 [-1 1]
% 输出参数：
% basis: 切比雪夫多项式基函数
function basis = chebyshev_basis(N, t)
if abs(t) > 1 + eps(max(1, abs(t))) * 10
    disp(t);
    error("chebyshev_basis_error t > 1");
end
    basis = zeros(length(t), N);
    for k = 1:N
        basis(:,k) = cos((k-1) * acos(t));
    end
end